Non-Abelian Topological Order in S'-Wave Superfluids of Ultracold Fermionic Atoms 
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It is proposed that in s-wave superfluids of cold fermionic atoms with laser-field-generated ef- 
fective spin-orbit interactions, a topological phase with gapless edge states and Majorana fermion 
quasiparticles obeying non-Abelian statistics is realized in the case with a large Zeeman magnetic 
field. Our scenario provides a promising approach to the realization of quantum computation based 
on the manipulation of non-Abelian anyons via an s-wave Feshbach resonance. 

PACS numbers: 



Introduction - Recently, there has been considerable 
interest in topological phases of quantum many-body 
systems, which are characterized by the following fea- 
tures P, [3, H i, i, H 0, 0: (i) there are topologically- 
protected gapless edge states on surface boundaries of 
the systems, which are stable against local perturbations, 
(ii) for two-dimensional (2D) systems, there are quasipar- 
ticles with fractional quantum numbers (e.g., fractional 
charges) termed "anyons" . To this time, the possibility 
of realizing topological phases has been studied for vari- 
ous states realized in condensed matter systems, such as 
quantum (spin) Hall states [2, K M, \ 4 , jja, M, vortex states 
of p + iv superconductors [^TbOUlll. , and spin liquid 
states I 7| . | l 3ll , and for cosmological systems such as axion 
strings [14j . The feature (ii) is particularly of interest in 
connection with the realization of fault-tolerant quantum 
computation based on the manipulation of non-Abelian 
anyons 0, E, EI EE 

UtI E8| . Since topological phases 
provide not only a novel paradigm of quantum ground 
states but also a potential breakthrough for technolog- 
ical advance, it is desirable to pursue various possible 
schemes for their realization. 

In this letter, we propose a scenario in which a topo- 
logical phase, possessing gapless edge states and non- 
Abelian anyons, is realized in a BCS s-wave superfluid 
(SF) of ultracold fermionic atoms in an optical lattice 
with a laser-field-generated effective spin-orbit (SO) in- 
teraction. It is possible to generate an artificial SO in- 
teraction that acts on atoms by using spatially varying 
laser fields El H EL "~ 



resonance [18j . However, unfortunately, this has not yet 
been realized because of huge loss [28] . Contrastingly, 
s-wave SFs of cold atoms with large BCS gaps have been 
realized via an s-wave Feshbach resonance |29l. l30| . Thus, 
our scenario based on s-wave SFs of cold atoms is deemed 
more advantageous for the realization of the topological 
order than that using p-wave SFs via a p-wave Feshbach 
resonance. Moreover, there is an important difference 
between the topological phase considered here and that 
of p-wave SFs. For chiral p-wave SFs, the non-Abelian 
anyons are vortices of the SF order parameter, which 
contain Majorana fermion modes. In striking contrast, 
in our system, the non-Abelian anyons are vortices of the 
SO interaction, i.e. the phase twist caused by the orbital 
motion accompanying spin flip. We propose an exper- 
imental scheme for generating and controlling vortices 
in the SO interaction, i.e. non-Abelian anyons, which 
are stabilized by use of a carefully designed laser setup 
rather than spontaneously formed macroscopic conden- 
sates. This scheme can be carried out by utilizing cur- 
rent sophisticated laser techniques. Thus, our proposal 
provides a promising approach to the realization of topo- 
logical quantum computation based on the manipulation 
of non-Abelian anyons. 

Model and Analysis of the topological phase - Let us 
consider an s-wave SF of neutral fermionic atoms in 
the 2D optical square lattice, which is described by the 
Hamiltonian H = 7Ykin + ^so + Us'- 



_ |22j, |2a| . The effective SO inter- 
action is a key factor in our scenario of the topologi- 
cal phase. Recently, topological phases in superconduc- 
tors and SFs have been investigated by several authors 
@, EE EE El EE US HE H3 . The previous studies, 
however, focus on the p-wave pairing state [Z(| • BCS gaps 
of p- wave superconductors in solid state systems are typi- 
cally very small, and it is difficult to utilize them for topo- 
logical quantum computation, because topological phases 
are destroyed by thermal excitations beyond bulk energy 
gaps. For cold atoms, in principle, a p-wave SF with a 

large BCS gap can be produced via a p-wave Feshbach where c| (ci a ) denotes a creation (annihilation) oper- 



ia p,=st,y ia 
i 

H SO = -^J2^4-xi c ^ -4+xi c n) 

i 

h s = -E^(4t c !i+ h - c -). 



2 



ator of the fermionic atom with pseudo-spin a = ("f , J.) 
at site i = (i x ,i y ), and -0s the gap function, x (y) is a 
basic lattice vector along the x(y)-axis. Tiso is an ef- 
fective Rashba type SO interaction [3l[. We will discuss 
later the method of generating the Rashba SO interaction 
for neutral atoms via laser fields. We also introduce the 
chemical potential fi and the Zeeman term induced by a 
magnetic field h. In the momentum space, the Hamilto- 
nian is recast into 
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,(3) 



where e& = —2t(cosk x + cosfc y ) — fi, gk = 
2A(sin/cj / , — sinfea;), and cr = (a x , <j y ) the Pauli matrices. 

As mentioned in the introduction, the non-Abelian 
topological order is characterized by the existence of gap- 
less chiral edge states propagating only in one direction 
and the existence of the non-Abelian anyons 0. The 
former is also associated with the nonzero Chern num- 
ber [33]. In the following, we demonstrate that these 
features are indeed realized in the system ([1]) when a cer- 
tain relation among /x, h, and ip s holds. (Eq.© below) 

A key observation of our analysis is that the Hamilto- 
nian Ti.(k) is unitary equivalent to the following "dual" 
Hamiltonian liP(k). 



ip s - ha z 
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with the unitary transformation 
H D (k) = DH(k)D\ D 
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(5) 



From Eq.([4|), it is found that the Rashba SO interac- 
tion gk ■ cr in the original Hamiltonian TL{k) is formally 
tranformed into a "p-wave SF gap" with the d vector, 
<ik = — gk, in the dual Hamiltonian HP(k). However, 
this does not necessarily mean that the topological prop- 
erties of Tt(k) are the same as those of a p-wave SF, 
since H D (k) has a non-standard constant kinetic term 
ej? = ip s . A similar p-wave SF state with a constant ki- 
netic energy term was considered before in the context of 
the quantum-Hall-effect (QHE) state Q. An important 
feature of |4]) is that the topological order emerges when 
/i, h, and ipg satisfy 



^ 2 + e(0,0) 2 </i 2 <^ 2 + e(^0) 5 



(6) 



with e(k x ,k y ) = tk- Here note that although the condi- 
tion ^ implies the Zeeman energy larger than the BCS 
gap ip s , the superfluidity is stable when A 3> h 33|, 34 1. 
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FIG. 1: The band energy of the lattice Hamiltonian |T} with 
edges at i x — and i x = 50(= L). Here k y G [— 7r,7r] denotes 
the momentum in the y-direction. We set t = 1, y, = —4, 
A = 0.5 and ip = 0.5. h is a) h = 0, b) h = 0.5, c) h = 0.8. The 
red thin line indicates a gapless chiral edge mode localized on 
the one side and green thick line a gapless chiral edge mode 
on the other side. They appear for \/tp 2 + e (0, 0) 2 < h < 

VV> 2 + e(o,^) 2 . 



This stability is specific to neutral atomic systems. For 
electron systems, such large magnetic fields usually de- 
stroy superconductivity via an orbital depairing effect. 

Let us first examine edge states in our model. Figure 
[T] illustrates the energy bands obtained by diagonalizing 
the lattice Hamiltonian (TT|) with the open boundaries at 
i x = 0, L for various h. Here we have taken the periodic 
boundary condition in the y direction, and k y £ [— tt, 7t] 
is the lattice momentum in the y direction. By increasing 
h from zero adiabatically, it is found that the bulk en- 
ergy gap closes at h — \/ips + e (0; 0) 2 (FigQJb)), then, 
for h satisfying ©, a gapless edge mode with a linear 
dispersion E ~ ck y [E ~ ~ck y ) localized on the one 
edge (the other edge) appears between the bulk energy 
gap (FigQJc)). This chiral edge state is stable against 
any weak local perturbations provided that there exists 
the nonzero Chern number; i.e. the topological num- 
ber equivalent to the total number of gapless chiral edge 
modes, which was first introduced in the case of the QHE 
states [32]. We calculated the Chern number Q for H(k) 
or equivalently 7i D (fc). (Since the Chern number is calcu- 
lated from the Berry curvature in the momentum space, 
it is not affect by the unitary transformation D which 
is independent of fc, ensuring the topological equivalence 
between (jSJ and j4}.) We found that Q = 1 when the 
condition ([6]) is satisfied [35] . This is consistent with the 
numerical results for edge states shown above. 

We, now, demonstrate that there exist the non-Abelian 
anyons in our system. For this purpose, we solve the 
Bogoliubov-de Gennes (BdG) equation for a single vor- 
tex: If there exists a single Majorana fermion zero mode 
for each vortex, vortices obey the non-Abelian statistics 
[3, We use the dual Hamiltonian 7i D to solve the 
BdG equation, then construct a solution in the original 
Hamiltonian TL by using the duality transformation ([5]). 
For simplicity, we assume e(0, 0) = for the time being. 
Then, low-energy properties are governed by fermions 
on the Fermi surface, which is split into |fe| ~ and 
|fe| ~ X/t by the SO interaction, but the larger Fermi 



3 



surface (|fc| ~ X/t) can be neglected for the zero mode 
[3^ . Thus, we concentrate on fermions with k w (0,0), 
for which 7i D (fe) is decomposed into the following two 
2x2 matrices, 



H D {k) = 



( tjj s — h 2X(k y + ik x ) 
\ 2X(k y - ik x ) -Vfc + h 
ip s + h 2X(—k y +ik x ) 
2X(—k y — ik x ) "V's — h 



(7) 



We consider a single vortex of the "p-wave SF gap" 
2X(±k y + ik x )- The BdG equations for Ti^_ with the sin- 
gle vortex can be solved by using the method developed 
in Then, we find a unique zero energy solution with 
a quasiparticle field 7^ = J dr[u ip\_ + v ip + ], where u = 

i ( re i6}-l/2 e -(h-4, s )r/2\ ^ Vo _ _ i r re -i6\-l/2 e -{h-il) i )r/2\ 

[35|] . The solution is normalizable when ([6]) is satisfied. 
This is the Majorana zero energy mode; i.e. 7T = 7. 
Using the duality transformation (O, we found that a 
vortex in the original Hamiltonian has a single Majorana 
zero mode given by D(uq, 0, vq, 0) t , which implies that 
the vortex is a non-Abelian anyon [9J . 

From the construction of the zero mode above, we no- 
tice that there is an important difference between a chiral 
p-wave SF and our system: While for a spinless chiral p- 
wave SF, a single Majorana zero mode exists in a vortex 
of the SF order parameter, for our non-centrosymmetric 
s-wave SF, a Majorana zero mode exists in a vortex twist- 
ing a phase of the SO interaction. This difference can 
be understood immediately from the duality (|4|) since a 
vortex in a gap function in the dual Hamiltonian is trans- 
formed into a vortex in the SO coupling. So far, we have 
assumed e(0, 0) = 0, but even when e(0, 0) 7^ 0, the ex- 
istence of the non-Abelian topological order is robust as 
long as h satisfies (6), because the topological character 
is not changed unless the bulk gap closes [25| ■ 

For the detection of the non-Abelian anyons, it is de- 
sirable that the zero energy state in a vortex is well sep- 
arated from excited states, the interaction with which 
may cause decoherence. The excitation energy in the 
vortex core is due to the kinetic energy which stems from 
the derivative term of the BdG equation, 2X( z fid y + d x ). 
Since the above solution for {uq,vq) indicates that the 
size of the vortex core is ~ 2X/(h — ip s ), the excitation 
energy is of the order ~ 2A/(2A/ (h—ip s )) ~ h—ip s . It can 
be tuned to be relatively large, and thus the experimental 
detection of the non-Abelian anyons is quite feasible. 

Possible Realization in Cold Fermionic Atoms - We 
now propose an experimental scheme for the realization 
of the topological phase mentioned above in ultracold 
fermionic atoms. It was recently pointed out by sev- 
eral authors that effective gauge fields interacting with 
atoms can be g enerated by spatially varying laser fields 
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z-direction [19]. The atoms occupy doubly-degenerate 
Zeeman levels of the hyperfine ground state manifolds, 
which are, respectively, the "spin up" state | |) and the 
"spin down" state | J,). We introduce the Zeeman field 
to lift the degeneracy. The Zeeman level split is denoted 
as Ez- It is assumed that standard tunneling of atoms 
between sites due to kinetic energy is suppressed by the 
large depth of the optical lattice potential. Tunneling 
of atoms between neighboring sites along the ^-direction 
(y = x, y) which conserves spins is caused by laser beams 
with the Rabi frequency via optical Raman transi- 
tions as proposed in Refs. 12, 36|; i.e. fl„o = Q'ii~l' 2 /2A 
with flj 2 the Rabi frequencies for the transition between 
the ground state and an excited state, and A detuning 
from the excited state. In addition, tunneling which ac- 
companies spin flip is also driven by two Raman lasers 
19]. In Fig[2l we show the optical lattice setup. The laser 
with the Rabi frequency fi„i (^2) is resonant for tran- 
sition I I) —> I J.) for the tunneling between neighboring 
sites in the forward (backward) ^-direction. As proposed 
in Ref . [l9| , the confining optical potential is tilted along 
both the ^-direction and the y-direction to assure that 
the forward and backward tunneling processes are respec- 
tively induced by the lasers with the different Rabi fre- 
quencies f2„i and Q^, which are required for the realiza- 
tion of the Rashba spin-orbit interaction as discussed be- 
low. The energy shift between nearest neighbor sites due 
to the tilting potential is A^ for the ir-direction and A y 
for the y-direction. We impose the condition A^, ^ A y to 
prevent tunneling with spin flip along the ?/(x)-direction 
due to the lasers with £l x (y)i,2- It is also assumed that the 
detuning from excited states for optical Raman transi- 
tions is much larger than A x ( y \, and thus the spatial vari- 
ation of the amplitudes of the Rabi frequencies due to the 
tilting potential is negligible. To realize the Rashba spin- 
orbit interaction for the two Zeeman levels, we choose the 
phases of the lasers as follows. The lasers are propagat- 
ing along the z-direction with an oscillating factor e lkzZ . 
The Rabi frequency Vt X 2 is expressed as Cl X 2 = \ fl x2 \e lkzZ . 
The phase of the laser fl x i is shifted by tt from that of 
£l X 2, and £l X 2 — — O x i holds. The phase of £l y \ (^2) is 
shifted by — tt/2 (tt/2), and Q y 2 — —iQ x ij Qy2 — — Q y i- 
Then, the laser-induced tunneling term which accompa- 
nies spin flip is expressed by Hso = J2Mx( c l- x i c i7 ~ 

C i+xl C n) + iX v( C L y l C i] - C \+yi C i]) + h - C -] With A " = 

c v J dripl(r - r l ^ fl )n„ 2 (r)i/j 1 (r - n), v = x, y, and 
c x = 1, c y = —i. Since we consider the 2D ccy-plane 
with z — 0, X v is real. For = X y , Hso is the Rashba 
SO interaction. 



These ideas can be utilized for our 
purpose. We consider fermionic atoms loaded in a 2D pe- 
riodic optical lattice, where there is no hopping along the 



To create vortices of the SO interaction, which are key 
ingredients for the realization of non-Abelian anyons, we 
replace the lasers that generate the SO interaction with 
those carrying orbital angular momentum parallel to the 
z-axis. Such lasers can be prepared by using Laguerre- 
Gaussian beams [371 . Then, a "vortex" of the SO interac- 
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FIG. 2: Setup of the confining optical potential for the ground 
state, v = x or y. Bold up and down arrows indicate, respec- 
tively, the spin-up and spin-down states. In this figure, ex- 
cited levels which mediate the hopping via two Raman lasers 
are not shown explicitly. 



tion is introduced: A„ — > \ v e %mB . The vorticity m is con- 
trolled by changing the configuration of the lasers. Fur- 
thermore, the spatially separated multiple vortices can 
be generated by using the following method. After in- 
troducing a vortex with vorticity m into the system, we 
switch the Laguerre-Gaussian beam to a Gaussian beam, 
i.e. a laser without angular momentum (37l |. The vortex 
still exists in the system, because of the conservation of 
the total angular momentum. However, for m > 1, the 
vortex with higher charges become energetically unstable 
toward dissociation into m vortices with single vorticity, 
and thus spatially separated multiple vortices are created 
in the system 33, 39| . This multiple vortex state allows 
the realization of the non-Abelian statistics of vortices. 

We can use the Feshbach resonance in the s-wave chan- 
nel for the formation of the s-wave Cooper pairs in this 
system [H, |3(|. Then, the topological phase described 
by the Hamiltonian (fT]) is realized. In this scheme, the 
role of the s-wave superfluidity is two-fold. One is to 
suppress bulk gapless quasiparticles which are harmful 
for topological stability. The other one is to generate the 
superposition of particles and holes, which results in the 
Majorana quasiparticles in vortices and edge states. It 
is noted that the vortices in the SO interaction are sta- 
bilized by the carefully designed laser setup rather than 
by macroscopic condensates. It should be emphasized 
that this experimental scheme is feasible for currently 
accessible laser techniques. As mentioned before, the 
non-Abelian anyons are stable for sufficiently low ener- 
gies <C mm{h — i/j a ,ip s }. Since ip s can be tuned to be 
large, i.e. ip s ~ Ep, by using the s-wave Feshbach reso- 
nance, the realization of the non-Abelian anyons in this 
scheme is quite promising. 

Summary- We have proposed a feasible scheme for the 
realization of the non-Abelian topological phase in an 
s-wave superfluid of cold atoms in an optical lattice, in 
which the non-Abelian anyons exist, and opens a possible 
way to realize topological quantum computation. 
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